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Scattering of hydrogen molecules from a reactive surface: 
Strong off-specular and rotationally inelastic diffraction 

Axel Gross and Matthias Schemer 
Fritz- Haber-Institut der Max- Planck- Ges ells chaft, Faradayweg 4-6, D-1^195 Berlin- Dahlem, Germany 

Six-dimensional quantum dynamical calculations of the scattering of H2 from a Pd(100) surface using a 
potential energy surface derived from density-functional theory calculations are presented. Due to the 
corrugation and anisotropy of the PES strong off-specular and rotationally inelastic diffraction is found. 
The dependence of the diffraction intensitities on the incident kinetic energy is closely examined. In 
particular we focus on the quantum oscillations for normal and off-normal incidence. 
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Coherent elastic scattering of atoms or molecules from 
surfaces is a tool for probing surface structures and 
adsorbate-substrate interaction potentials. In particu- 
lar helium atom scattering (HAS) has been used inten- 
sively to study surface crystallography and the shape 
of physisorption potentials (see, e.g., pj| and references 
therein). Hydrogen molecules has been utilized less fre- 
quently in order to study interaction potentials. The 
coherent elastic scattering of molecules is more complex 
than atom scattering because in addition to parallel mo- 
mentum transfer the internal states of the molecule can 
be changed in scattering. This means in particular ro- 
tational transitions, because the energies typically em- 
ployed in molecular scattering experiments are not suffi- 
cient to excite molecular vibrations. 

Rotational transitions in elastic coherent scattering 
lead to the appearance of additional peaks in the diffrac- 
tion pattern since the total kinetic energy of the molecule 
is changed. Most experimental studies of rotationally in- 
elastic scattering of hydrogen have been performed at the 
surfaces of ionic solids ||[|. (Here "inelastic" refers to 
energy exchange between the hydrogen degrees of free- 
dom, not with respect to the substrate.) In the hitherto 
existing experiments of hydrogen scattering from metal 
surfaces the additional peaks could usually hardly be re- 
solved except in the case of HD scattering ||f7|,[| , 
where the displacement of the center of mass from the 
center of the charge distribution leads to a strong rota- 
tional anisotropy. Only recently Bcrtino et al. [|| clearly 
resolved rotationally inelastic peaks in the diffraction 
pattern of D2/Ni(110) in addition to rotationally elastic 
peaks, which had been observed before by Robota et al. 
||. Furthermore, Cvetko et al. reported diffractive 
and rotationally inelastic scattering of D2/Rh(110). 

In the case of hydrogen scattering from reactive sur- 
faces, the repulsive interaction is not mediated by the 
tail of the metal electron density, but occurs rather 
close to the surface Due to the chemical nature 

of this interaction the potential is strongly corrugated 
and anisotropic with regard to the molecular orienta- 
tion. This requires a realistic microscopic description of 
the scattering dynamics. Existing theoretical studies of 



hydrogen scattering at surfaces have either treated non- 
reactive systems [ p^| , p^l5|jr^ , p7[ | or they were restricted 
to low dimensions on model potentials |18 19| , |20[ |. 

We have performed six-dimensional quantum dynam- 
ical calculations of the scattering of H2 from Pd/(100) 
on a PES that was derived from density-functional the- 
ory calculations pd| , ^l| . The PES exhibits a coexis- 
tence of activated with non-activated paths to adsorp- 
tion together with a broad distribution of barrier heights. 
All six degrees of freedom of the hydrogen molecule 
are treated quantum dynamically in a coupled-channel 
scheme |^2|,^3| while the substrate atoms are kept fixed; 
for computational details we refer to Ref. |2^| . 

Our results show strong off-specular and rotationally 
inelastic diffraction which are caused by the corrugation 
and anisotropy of the PES. The intensities of the diffrac- 
tion peaks exhibit a pronounced oscillatory structure as 
a function of the incident energy. We have already pre- 
dicted such an structure for the sticking probability of 
H 2 /Pd(100) [^T|. Very recently Rettner and Auerbach 
searched for these oscillations experimentally [p5| , but 
did not find any. We will carefully discuss the physical 
origin of the oscillations and point out why they might 
not have been detected. Hence this paper may also be 
regarded as an extended comment to Ref. p^| . 

One typical calculated angular distribution of H2 
molecules scattered at Pd(100) is shown in Fig. [j]. The 
total initial kinetic energy is Ei = 76 meV. The inci- 
dent parallel momentum corresponds to 2hG along the 
(Oil) direction, where G = 2ir/a, a is the nearest neigh- 
bor distance between two Pd atoms. This leads to an 
incident angle of 9i — 32°. The molecules are initially 
in the rotational ground state ji — 0. (m, n) denotes 
the parallel momentum transfer AG|| = (mG,nG). The 
specular peak is the most pronounced one, but the first 
order diffraction peak (10) is only a factor of four smaller 
(note, that in a typical HAS experiment the off-specular 
peaks are about two orders of magnitude smaller than the 
specular peak). The results for the rotationally inelastic 
diffraction peaks j = — > 2 have been summed over all 
final azimuthal quantum numbers mj, as has been done 
for all results presented in the following. The excita- 
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FIG. 1. Angular distribution of scattered H2 molecules. 
The initial kinetic energy is Ei — 76 meV, the incident angle 
is 8i = 32° along the (Oil) direction. The molecules are ini- 
tially in the rotational ground state ji — 0. The rotationally 
inelastic diffraction peaks have been marked with Aj = 2. 



tion probability of the so-called cartwheel rotation with 
rn = is for all peaks approximately one order of mag- 
nitude larger than for the so-called helicopter rotation 
to = j, since the polar anisotropy of the PES is stronger 
than the azimuthal one. 

This steric effect in scattering could already be ex- 
pected from detailed balance arguments. According to 
the principle of detailed balance, in a equilibrium situa- 
tion the flux impinging on a surface from the gas-phase, 
which is rotationally isotropically distributed, should 
equal the desorption plus the scattering flux. Since in 
desorption the helicopter rotations are preferentially oc- 
cupied |^l],^6| , in scattering the cartwheel rotations have 
to be preferentially excited. 

The intensity of the rotationally inelastic diffraction 
peaks in Fig. [j] is comparable to the rotationally elas- 
tic ones. Except for the specular peak they are even 
larger than the corresponding rotationally elastic diffrac- 
tion peak with the same momentum transfer (to, n). Note 
that due to the initial conditions the rotationally elastic 
and inelastic (20) diffraction peaks fall upon each other. 

The out-of-plane scattering intensities are not negligi- 
ble, which is demonstrated in Fig. 0. The initial condi- 
tions are the same as in Fig. |. The open circles cor- 
respond to rotationally elastic, the filled circles to rota- 
tionally inelastic diffraction. The radii of the circles are 
proportional to the logarithm of the scattering intensity. 
The sum of all out-of-plane scattering intensities is ap- 
proximately equal to the sum of all in-plane scattering 
intensities. Interestingly, some diffraction peaks with a 
large parallel momentum transfer still show substantial 
intensities. This phenomenon is well known form helium 
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FIG. 2. Angular distribution of the in-plane and 
out-of-plane of scattering of H2/Pd(100). The initial con- 
ditions are the same as in Fig. hi. Open circles correspond 
to rotationally elastic, filled circles to rotationally inelastic 
diffraction. The radii of the circles are proportional to the 
logarithm of the scattering intensity, x denotes the (Oil) di- 
rection, y the (011) direction. The specular peak is the largest 
open circle. 



atom scattering and has been discussed within the con- 
cept of so-called rainbow scattering p7fl . 

A more detailed analysis of the diffraction intensities 
versus total incident kinetic energy is shown in Figs. ^M. 
In Fig. H the intensities of four different diffraction peaks 
are plotted for normal incidence; Fig. ||a) shows rotation- 
ally elastic, Fig. ||b) rotationally inelastic diffraction. In 
the experiment the molecular beams are not monoener- 
getic but have a certain velocity spread. In order to make 
close contact to experiments we assumed for our results 
in Figs.|||-||an initial velocity spread Avi/vi = 0.05, close 
to the velocity spread of recent experiments || . 

The theoretical data still exhibit a rather strong oscil- 
latory structure which is a consequence of the quantum 
nature of hydrogen scattering. We first focus on the spec- 
ular peak in Fig. |3^i). The oscillations at low kinetic en- 
ergies reflect the opening up of new scattering channels 
p8[ |. The first pronounced dip at Ei — 12 meV coin- 
cides with the opening up of the (11) diffraction peak, 
the small dip at Ei = 22 meV with the opening up of the 
(20) diffraction channel. The huge dip at approximately 
50 meV reflects the threshold for rotationally inelastic 
scattering. Interestingly, the rotational elastic (10) and 
(11) diffraction intensities show pronounced maxima at 
this energy. This indicates a strong coupling between 
parallel motion and rotational excitation. 

At larger energies the scattering intensities still show 
an oscillatory structure. In addition to the opening up 
of new scattering channels this also reflects the existence 
of scattering resonances E9] : molecules can become tem- 
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FIG. 3. Scattering intensity versus kinetic energy with an 
initial velocity spread Avi/vi — 0.05 for H2 molecules in the 
rotational ground state under normal incidence, a) Rotation- 
ally elastic j = — > scattering, b) rotationally inelastic 
j = —* 2 scattering. The intensities are summed over all 
final azimuthal quantum number rrij. 



porarily trapped into metastable molecular adsorption 
states at the surface due to the transfer of normal mo- 
mentum to parallel and angular momentum which reso- 
nantly enhances the scattering intensity. 

These structures are known for a long time in He 
and H 2 scattering (2| and also in LEED For 
H 2 /Pd(100), however, measuring these oscillations is a 
very demanding task. They are very sensitive to surface 
imperfections like adatoms or steps. Due to the high reac- 
tivity of this system either an adlayer of hydrogen atoms 
builds up very rapidly during the experiment or in order 
to avoid this one has to go to high surface temperatures 
|p2| . Both suppresses the oscillations. 

Figure |^b) shows the intensities of rotationally inelas- 
tic diffraction peaks. The specular peak is still the largest 
one, however, some off-specular peaks become larger than 
the (00) peak at higher energies. This is due to the fact 
that the rotationally anisotropic component of the poten- 
tial is more corrugated than the isotropic component jll] . 
Furthermore it is apparent that the oscillatory structure 
is somewhat smaller for rotationally inelastic than for ro- 
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FIG. 4. Scattering intensity versus total kinetic energy 
with an initial velocity spread Avt/vi = 0.05 for H2 molecules 
in the rotational ground state under off-normal incidence. 
The incident parallel momentum corresponds to 2hG along 
the (Oil) direction with G — 2n/a. The incident angle is 
given at the top axis, a) Rotationally elastic j = — ► 
scattering, b) rotationally inelastic j — — > 2 scattering. 
The intensities are summed over all final azimuthal quantum 
number m, . 



tationally elastic scattering. 

The results shown in Fig. [| correspond to a H2 beam 
impinging on the Pd(100) surface under off-normal in- 
cindence. The initial parallel momentum is pii = 2hG 
along the (011) direction. Note that the incident angle 
decreases with increasing total kinetic energy for con- 
stant initial parallel momentum. Only in-plane scatter- 
ing intensities are plotted. Now the intensities of the 
off-specular peaks compared to the specular peak are 
much larger than in the case of normal incidence, which 
means that the change of parallel momentum in molecule- 
surface scattering is much more probable for off-normal 
than for normal incidence. Again a pronounced oscilla- 
tory structure for the single diffraction peaks is observed. 

One overall trend is evident in Fig. ^ and ^: the scat- 
tering intensities of almost all single peaks decrease with 
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FIG. 5. Total scattering intensity versus total kinetic en- 
ergy with an initial velocity spread Avi/vi = 0.05 for ro- 
tationally elastic and inelastic scattering at normal inci- 
dence and off-normal incidence with initial parallel momen- 
tum p|i = 2hG. 



increasing kinetic energy as also observed in the exper- 
iment ||, whereby the specular peak decreases more 
strongly than the off-specular peaks. In ref. JoJ the ener- 
getic dependence of the scattering intensities is discussed 
within a Debye- Waller analysis in the impulsive and har- 
monic approximation jlj ; the decrease in the peak height 
with increasing energy is solely attributed to the inter- 
action with surface phonons. Our analysis, like others 
before J^,|l9|,^0|] , shows that there is already a substan- 
tial decrease in the intensities for hydrogen scattering at 
metal surfaces without taking any energy transfer to the 
phonons into account. 

In Fig. U we have plotted the total reflection prob- 
ability for rotationally elastic and rotationally inelastic 
j = — ► 2 scattering, i.e. the sum over all corresponding 
diffraction peaks. After an initial increase, which corre- 
sponds to the decrease of the sticking probability at low 
energies [^l],[53| , the reflection probabilities generally de- 
crease with increasing energy. For kinetic energies above 
0.1 eV the probability of rotationally inelastic j = — > 2 
scattering becomes larger than the probability of rota- 
tionally elastic scattering. In addition, the rotational ex- 
citation is larger for off-normal incidence than for normal 
incidence. The oscillatory structure of the total scatter- 
ing intensities is less pronounced (note the different in- 
tensity scale) than for the single diffraction peaks, but is 
still clearly visible, especially the threshold for rotational 
excitation at normal incidence. 

Note that the minimum in the specular peak at 
Ei w 12 meV has now turned into a local maximum in 
the total scattering intensity. That is due to the fact that 
at this energy not only one scattering channels opens up, 



but actually four due to symmetry, namely the (10), (10), 
(01) and (01) peaks. The additional scattering intensity 
in these diffraction peaks over-compensates for the loss 
of the specular peak in the total scattering intensities. 

Although the scattering intensities of the single peak 
for off-normal incidence also showed a strong oscillatory 
structure, their sum is rather smooth (see Fig. ||). This 
can also be understood by symmetry arguments. For an 
incident angle along the (011) direction the (10), (10), 
(01) and (01) peaks, e.g., open up at three different ener- 
gies, for a general incident azimuth even at four different 
energies. Therefore the effect of the opening up of new 
scattering channels is significantly reduced for the total 
scattering intensity compared to normal incidence. As 
far as quantum oscillations are concerned, it is therefore 
not reasonable to compare experiments for an angle of 
incidence of 9i = 15° [^5| with calculations for normal 
incidence. 

In order to observe quantum oscillations we suggest 
not to determine the whole scattering flux for off-normal 
incidence as done in ref. |25f| , but rather monitor the in- 
tensity of a single diffraction peak as a function of the 
kinetic energy. Thus one also avoids to include inco- 
herently scattered molecules due to impurity or multi- 
phonon scattering in the measurements. In our calcula- 
tions the substrate is held fixed. Substrate motion cer- 
tainly makes the detection of quantum oscillations more 
difficult. Since these oscillations have been observed for 
H 2 scattering at other surfaces p9| , they should, how- 
ever, also be detectable for H 2 /Pd(100) scattering. 

The general decrease of the total scattering intensities 
at higher energies reflects the fact that the sticking prob- 
ability rises in this energy range so that less molecules 
are reflected. However, the decrease of the total intensi- 
ties is much slower than the decrease of the single peaks 
in Figs. H and ^ (note again the different intensity scale). 
This implies that the decrease in the intensities of the 
single peaks is partly compensated by the opening up of 
new scattering channels with increasing kinetic energy. 
The number of energetically accessible parallel momen- 
tum states increases roughly linearly with energy, as does 
the number of rotational states, so that the number of ac- 
cessible scatterings states increases quadratically with in- 
creasing energy (again neglecting vibrational excitation) . 
Hence the increase of the available phase space for scat- 
tering with increasing kinetic energy is responsible for 
the stronger decrease of the single diffraction peaks com- 
pared to the total reflection probability. 

In conclusion, we have performed a six-dimensional 
quantum dynamical study of the scattering of H 2 from a 
Pd (100) surface using a potential energy surface (PES) 
derived from first-principles calculation. The corrugation 
and anisotropy of the PES lead to strong off-specular and 
rotationally inelastic diffraction intensities. Strong quan- 
tum oscillations in the diffraction intensities as a function 
of the kinetic energy are obtained. Especially at low ki- 
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netic energies they reflect the opening up of new scatter- 
ing channels with increasing kinetic energy. 
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